Is the QCD ghost dressing function finite at zero momentum 
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Qh' Abstract 

< 

^O ■ We show that a finite non-vanishing ghost dressing function at zero momentum satis- 
fies the scaling properties of the ghost propagator Schwinger-Dyson equation. This kind 
of Schwinger-Dyson solutions may well agree with lattice data and provides an interest- 

\j^ ' ing alternative to the widely spread claim that the gluon dressing function behaves like 

l/~j . the inverse squared ghost dressing function, a claim which is at odds with lattice data. 

^D ' We demonstrate that, if the ghost dressing function is less singular than any power of 

'^ . p, it must be finite non-vanishing at zero momentum: any logarithmic behaviour is for 

^^ ' instance excluded. We add some remarks about coupled Schwinger-Dyson analyses. 
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^ ■ 1 Introduction 



The infrared behaviour of Landau gauge lattice gluon and ghost propagators is an interesting 
and hot subject. Two main methods are used: lattice QCD (LQCD) and Schwinger-Dyson 
equations (SDE) in which we include related methods as RGE, etc. In ref. Q we have shown 
that a combination of both methods is extremely enlightening as it combines the advantages 
of lattice QCD's full control of errors and SDE's analytical character. 
We only consider the particularly simple ghost propagator SDE : 
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We have studied the discrepancy between LQCD data and a widely spread beUef: the 
ghost propagator SDE is claimed to imply a gluon dressing function behaving like the inverse 
squared ghost one. In ref. ^ we have reconsidered the scaling properties of the SDE and 
found three possible ways out of this problem (which are summarised in table 1 of that paper). 
The first one is to assume a singular behaviour of the ghost-gluon vertex in the deep infrared. 
A second possibility implies a very singular ghost dressing function which is excluded by 
LQCD. The third one is to assume that the ghost dressing function is less singular in the 
infrared that any power of p. In view of the general belief that the ghost dressing function 
was strongly singular we had not paid in ref. ^ attention to the third one. 

Very recently, Sternbeck et al. |2j have produced two new evidences: i) the ghost-gluon 
vertex seems not to be singular, ii) the ghost dressing function seems to behave at most like 
logp in the infrared. These two evidences, taken together, strongly encourage us to consider 
now seriously the third above-mentioned solution. This is the aim of the present letter. 

To our surprise we found that one can demonstrate from the scaling analysis of the ghost 
propagator SDE the impossibility of a logp behaviour or any other behaviour which is less 
divergent than any power of p: under these conditions, the ghost dressing function necessarily 
has a finite non-vanishing limit at zero momentum. This is at odds with a very general belief 
that the ghost dressing function is divergent. The proof will be displayed in section |31 We 
will shortly discuss published results about coupled gluon and ghost SDE in sectional 

2 Notations and summary of up-to-date lattice results 

We use the following notations p^: 

ff,{-q,k;q-k) = q^Hi{q,k) + {q - k)^H2{q,k) 

where G^"^' and F'^'^' are respectively the gluon and ghost propagators, G and F are respec- 
tively the gluon and ghost dressing functions and where V^{—q,k]q — k) is the ghost-gluon 
vertex ( k and q are the momenta of the incoming and outgoing ghosts and q — k the mo- 
mentum of the gluon) . 

Following for simplicity the common, convenient, but not really justified, assumption of 
a power-law behaviour of the propagators in the deep infrared, we define 

F{k^) ^ (^-^ ,G(fc2)~(^j . (2) 

In ref. ^ we have also defined an infrared exponent a-p for the vertex funtion Hi (ar < 
means a singular infrared behaviour). 

Using the ghost propagator SDE equation it is often claimed that 2aF + ac = 0. This 
belief is so strong that one often uses only one parameter k = —ap = Q-gI'^- However, as 
we will see in more details, everybody agrees that ac is close '^ to 1 and it becomes now 
clear ^ that otp is close to zero. Then the relation 2ap + ac = is not satisfied \^ and 
the arguments which support it have to be reconsidered. 

■^For example, in many SDE approaches it is found Pj ao — 1-18 

*One may wonder why many power law fits have given negative of- Our own fit in ref.0 (table 2) has 
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Figure 1: F{p) from lattice simulations for SU{2) (left, Vol = 48^, P = 2.3) and SU{3) 
{rightjV ol = 32^, /3 = 5.8). /3 = 2.3 for SU{2) has been chosen to guarantee that the string 
tension in lattice units is close to that of /? = 5.8 for SU{3). 

The lattice gluon propagator. Several SDE studies ([3] and references therein) predict 
a vanishing zero momentum propagator while, as discussed in ref. [1], a gluon propagator 
converging continuously to a non-zero value at vanishing momentum is a rather general lattice 
result (in particular, the authors of ref. |H1 obtain a non- vanishing infrared limit for the gluon 
propagator at a lattice volume of around 2000 fm^). Therefore our preferred solution ^ is 

ac = 1. (3) 

But, even if we relax this relation and assume a vanishing gluon propagator with ac > 1, 
the solution with a finite ghost dressing function at zero momentum remains possible as we 
shall see. 

The lattice ghost propagator. Very recent lattice estimates [2] seem to point towards 
a ghost dressing function rather close to the perturbative behaviour: the dressing function 
only shows, if any, a logarithmic dependence on the momentum (see fig. 2 of ref. [2|). 

We confirm these results. In Fig. ^ the ghost dressing function is plotted as a function of 
log{p) for small values of the momenta. These plots were obtained from lattice simulations 
at /5 = 5.8 and a volume 32^ in the SU{3) case and at /3 = 2.3 and a volume 48^ in the SU(2) 
case. It is clear from these plots that F{p) does not exhibit any power law: ap = 0. For 
SU(3) F{p) is approximately linear in log(p) and for SU{2) it has even a smoother behavior 
(In this case one obtains a good fit of the data with a function C(log \p\y' and 7 ~ 0.4). 



The ghost-gluon vertex. In ref. J2], the authors did not find any evidence for a singularity 
in the case of a vanishing gluon momentum. Let us remark that this particular kinematical 
configuration isolates the form factor Hi (see Eq.(^) which enters in the SD equation (It 

produced negative values very close to zero, but the errors were clearly underestimated. Presumably the 
systematic one, due to the functional form chosen for the fit, has not been properly taken into account. This 
is also the case in several other published results. 

^Let us recall however that there is still a problem coming from the Slavnov- Taylor identity for the three 
gluon- vertex: we have shown that, the vertices being regular when one momentum tends to zero, it implies 
aa < 1- Is it possible to avoid any contradiction by assuming, as done by Cornwall [11| . a non-regular 
behaviour for the longitudinal part of the three-gluon vertex? This deserves more investigation. 



is worth recalling that in perturbation theory Hi{q,0) + H2{q,0) is equal to 1 although 
Hi{q,0) is not [Jj)- Our dimensional analysis of the ghost SDE (see next section |2I) invokes 
a different kineniatical configuration for this form factor, Hi{q,k -^ 0) instead of Hi{k,k). 
The non-singular behaviour they found as k tends to excludes however the singularity of 
the ghost-gluon vertex that we proposed in ref. [l] as our favoured solution to reconcile the 
ghost propagator SDE and the lattice inspired relation 2aF + ac > 0. 

In conclusion lattice simulations show a strong evidence that ac = 1 and 2aF + ac > 0, 
far from zero. Now we have a fair indication that ap = and about the regularity of the 
vertex form factor involved in the ghost SDE. This leads us to revisit, in next section, the 
case in column 4 of table 1 in ref. ^ for ghost and gluon propagators and vertices satisfying 
the scaling properties of ghost SDE. 

3 Ghost SDE: the case ap = 

We will now demonstrate that F{0) is finite non-vanishing for ap = 0. We will exploit the 
constraints, summarized in table 1 of ref. P, between the infrared exponents ap, ac, ctr, 
from the ghost SDE. The IR convergence of the loop integral in the ghost SDE implies the 
two conditions: 

ap + ar > —2 , og + ar > —I ■ (4) 

Then the dimensional consistency of ghost SDE at small momenta leads to only three allowed 
cases: 

i) ap 7^ and ap + ac + ar < 1 =^ 2ap + ac -|- ar = 

ii) ap 7^ and ap -\- aQ -|- ar > 1 ==^ ap = — 1 

iii) ap = and ac + ar ^ 1 does not require any further constraint. 

We shall look, in the following, at the consequences of the third case. It includes in 
particular ac = 1 and ar = which is favoured by lattice simulations (see section ^ ). 
Nevertheless we shall not suppose, in the following derivation, anything more than ag+ar > 1 
and conditions Q. This leaves open, for example, the possibility that the gluon propagator 
goes to zero in the IR limit, the vertex remaining finite or singular. 

Of course, even with ap = 0, we cannot exclude a priori the possibility that F[k) diverges 
or tends to zero more slowly than any power of k when /c — > 0. We shall however prove that 
this is not allowed: F{k) remains finite in this limit provided that the two following conditions 
are satisfied: 



ap = , ac + ar > 1 (5) 

Writing the subtracted bare SD equation for two scales \k and nXk (see Eq.(ll) of ref. P) 
one obtains: 



1 ,i^[J^(E^(^±^_^^ 
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F{Xk) F{KXk) "^ V (27r)4V g2 \^ ^2 
-G{{q - \kf)Hl{q, \k) G{{q - ti\kf)Hl{q, nXk) 



{q - Xky {q - KXkf 



(6) 



where A is a parameter which we shall use to study the IR (A -^ 0) dimensional behaviour 
of -F; K is a fixed number, < k < 1, needed to write a subtracted equation ensuring the 



UV convergence. It was shown in ref. ^ that the r.h.s. of Eq. Q is the sum of two terms 
behaving respectively as ^2Mm{oF+OG+"r,i) ^jj^j ^2 -^j^gn _\ _!. q, go it behaves as A^ when 
the conditions © are satisfied. For any k there is a value of A and c such that VA' < A we 
liave l^p^ - :p^| < cA'2, thus: 



1 1 



'F{Xk) F{nXk) 
1 1 



< cX' 



'F{K^-^\k) F{K^\k) 
which implies: 

1 1 , 1 - K' 



< cA2k2("-i) (7) 



2n 

< c-^X\ (8) 



'F(A/c) F{K-^Xk)' - 1-K 

So -F — > 00 when A — > is excluded because taking the limit of the above expression when 
n — > 00 we should have Ipvi^l < Ct^-^A^ and F would diverge as or more rapidly than p- 
implying ap < — 1 in contradiction with the hypothesis ap = 0. Let us remark that F ^ 
is also excluded: Eq. (jH)) implies I piJ^Xk) I — I F(\k) I "*" c \"l'^a A^ and p/Jaxk) cannot tend 
to infinity when n —> 00. This completes the proof. Notice that we have used bare Green 
functions and couplings, everything remains however exactly the same if we replace them by 
renormalized ones. 

This is our main result: If ap = the ghost dressing function has to be finite 
and 7^ in the IR limit. This solution is compatible with our knowledge from lattice 
simulations about the behavior of the ghost dressing function and ghost-gluon vertex. Of 
course, the current lattice simulations cannot yet exclude a smooth divergence which the 
preceding dimensional analysis forbids. A detailed numerical study of the ghost propagator 
in the deep IR is strongly needed. 

4 Remarks about coupled gluon and ghost SDE solutions 

The combination of the scaling analysis of ghost SDE and lattice predictions appears to be 
very restrictive concerning the low-momentum behaviour of gluon and ghost propagators. 
Such a behaviour must be a solution of the combined SDE for both gluon and ghost propaga- 
tors. The schemes followed to solve the combined SDE's have often led to 2ap + ac = and 
ac ^ 1 (hence a strongly divergent ghost dressing function). However, a two- loop analysis [H] 
proves to be much less restrictive in constraining ac- Our findings require to reconsider these 
approaches by taking into due account the special case ap = 0. 

In a recent paper ^01 Aguilar and Natale found ac = 1 and ap = —0.04, not far from 
our present conclusions and deserving a closer comparison. They followed the Cornwall ^J 
prescription for the trilinear gluon vertex and solved the coupled equations for the ghost and 
gluon propagators in the Mandelstam approximation J21- Concerning the ghost dressing 
function, in spite of the fact that they find it slowly power-like divergent, it remains flat till 
very small momenta. This last point is in contradiction with the lattice results in ref. |2] 
and ours in Fig. ^ where F{k) is not at all so flat and shows a logarithmic enhancement as 
the momentum decreases. Of course, if power-like divergences are excluded, the arguments 
presented in section |21 imply a flat dressing function in a small momentum range presumably 
not yet reached by current lattice analyses. 
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Table 1: best- fit parameters. 

To compare quantitatively their gluon propagator with LQCD 0, we have applied the 
simple parametrization they proposed: 
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q^ + mo{Ly^ 



where a stands for the lattice spacing and L for the lattice length. In ref. [2], the gluon 
propagator was estimated from 24^ lattices at /? = 5.6, 5.8, 6.0 and 32^ lattices at /? = 5.7 and 
P = 6.0 and analyzed through an instanton liquid model that failed in describing the very 
low momentum range {q < 0.4GeV). In fig. [21 we plot the curves corresponding to the best-fit 
parameters niQ and Z(, collected in tabled The parametrization Eq. @ matchs pretty well 
the lattice data ^. Moreover one knows that, at the leading log. 
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d{log{Zb{a)) = -d{\ogP) . 



(10) 



Performing a linear fit of log(Z{,) as a function of log(/3) for /? > 5.7 we obtain a slope approx- 
imately equal to 0.69 which has to be compared to gf = 0.59. That result is unexpectedly 
good for the large lattice spacings we take in consideration. 
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Figure 2: Best fits to the lattice data (left) and Log-log plot of Zi, in terms of the lattice bare gauge 
coupling parameter (3 (right). The solid blue (dotted red) line shows a fit to a linear formula where 
the slope is free to be fitted (fixed by one- loop perturbation theory in Eq. (|10|l ). 



^The masses we obtain differ from the one quoted in ref. |1()) but these depend on a parameter, A, which 
in their approach can be varied. 



5 Conclusions 

The main result we presented in this brief note was to emphasize the interest of a general class 
of SDE solutions where ap ~ 0. This solution has the advantage of being compatible with 
other convincing lattice results namely: ac = 1 and 2aF + aG > 0. It is also compatible with 
a still uncertain result concerning the regularity of the ghost-gluon vertex. We have proven 
that if ap = the ghost dressing function must be finite in the IR limit. Of course one would 
need measures on larger volumes in order to test the finiteness of the ghost propagator in the 
limit k ^ 0. 

We have discussed some results from published coupled ghost and gluon SDE solutions 
and also shown that the lattice gluon propagator data at low momenta can be well described 
by the very simple parametrisation Eq. @ inspired by a recent gluon SDE analysis. 
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